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The Volatility in a Multi-share Financial Market Model 
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Abstract. Single index financial market models cannot account for the empirically observed complex in- 
teractions between shares in a market. We describe a multi-share financial market model and compare 
characteristics of the volatility, that is the standard deviation of the price fluctuations, with empirical 
characteristics. In particular we find its probability distribution is similar to a log normal distribution but 
with a long power-law tail for the large fluctuations, and that the time development shows superdiffusion. 
Both these results are in good quantitative agreement with observations. 

PACS. 5.40 Fb Random walks and Levy flights - 87.23 Ge Dynamics of social systems 



1 Introduction 

Recent works have shown that financial markets cannot be 
completely described by single index models since they do 
not account for complex interactions among stocks. Share 
and currency cross correlation matrices show some large 
eigenvalues where the market and certain groups of com- 
panies/currencies move together, against a background 
' that would be expected from Random Matrix Theory B 
, The largest eigenvalue shows sudden increase at crashes 
showing a strong global behaviour at such times. Share 
cross-correlations are generally positive jj| . Recent results |i 
show more complex properties of the simultaneous distri- 
bution of individual share returns coupled in a market. 
These results indicate that individual stocks in a market 
are coupled in a complex fashion and single index models 
cannot account for such characteristics. 

In a recent paper || we have described a multi-share 
model of a financial market and shown that it can ac- 
count for some inter-share characteristics as well as for 
some of the now 'well-known' properties, such as the ob- 
servation that real market returns distributions show 'fat- 
tails'. That is, the returns distributions are Levy stable 
distributions for the central part often characterised by 
a parameter a ~ 1.4 — 1.7, while the wings (after about 4 
standard deviations) fall-off faster as a power-law (expo- 
nent 2^5) or possibly as a streatched exponential 

In this paper we will show by numerical simulations 
that this model can also account for the empirically ob- 
served properties of the volatility. Volatility is the local 
standard deviation of a price returns time series. The 
volatility distribution shows an approximate log-normal 
distribution for the central part but with tails that seem 
to follow a power law with exponent ~ 4,|Rq|. Further- 



more volatility is known to be 'clustered'. That is we 
observe bursts of relatively high volatility separated by 
longer periods of relatively low volatility. This intermit- 
tent behaviour is characterised by a superdiffusion with 
exponent ~ 0.7»[p|,|TT 



2 Model 

Here we describe the model first presented in ju| . A more 
detailed discussion of its origins can be found in || . Pre- 
cursor models can be found in @H,|||. 

In this model there are N interacting stocks i = 1, N 
where for example i — 1 refers to 'IBM', i — 2 refers 
to 'Disney' etc. If we were considering the S&P500 then 
N » 500. Each stock is completely characterized by two 
variables: 

1) Excess demand/supply, Si(t). This is a spin vari- 
able Sj(i) = ±1 which describes stock i's current market 
demand/supply state. That is if Sj(t) = 1 then stock i is 
in the excess demand state at time t, but if Si(t) = — 1 
then it is in the excess supply state. In reality excess de- 
mand/supply has magnitude as well as sign but in this 
extreme simplification only its sign is taken into account. 
This is because as explained below we are only interested 
in modeling the self-reinforcing persistency characteristics 
of each stock i. (This model can be easily generalized by in- 
cluding a excess supply/demand magnitude gi(t), by, say, 
randomly choosing the magnitude gi{t) from a Gaussian 
(or otherwise) distribution each t.) The price return Api(t) 
for stock i at time t is given by Api(t) — Nsi(t)/2. The 
'index' price return Ap{t) is then given by Ap(t) = S(t)/2 

where S(t) = X)i=i s i(*)- This is just the usual linear re- 
lationship dp/alt = Dem{t) — Sup(t), where Dem(t) and 
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Supit) refer to the demand and supply at time t respec- 
tively. (See for example fl8|| ). 

2) 'Speculator Confidence', Viit) in stock i at time t. 
This is a real number which describes the overall (mean) 
speculator expectation of the current excess demand/supply 
state Si(t) continuing to the next time step, or reversing. 
In this model speculators are homogenous in their views. 
High confidence in a stock means speculators in general 
expect the current demand/supply state Sj(i) of that stock 
to continue, low confidence means they expect a reversal. 
Our confidence Vi (t) is cumulative and should be thought 
of like a 'fitness' in ecodynamics. It is not an intrinsic 
property of a stock, but is defined only by the current 
speculators, as a result of their expectations, and of course 
changes with time. 

There is no spatial dimension and stocks interact only 
through the market macrostate defined by 1) the market 
excess supply/demand G(t) = S(t)/N where G(t) is the 
mean-spin G(t) = l/JV^ti (t), and 2) the market con- 
fidence (market fitness) defined by V(t) = 1/N YliLi Vi CO- 
In real markets both of these are known (imperfectly) to 
the speculators. (In the case of the mean-confidence V(t), 
there may just be a general 'mood' traders can sense.) 

For the dynamics every time step t we calculate G(t) 
and V(t) and then calculate the relative confidence Ui(t) of 
stock i given by Ui(t) = Vj(t) — V{t). In this model specula- 
tors' expectations are on average completely self-fulfilling 
and with probability Qi(t) = Q{ui{f)) where Q(x) is given 

by, 



Q(x) 



1 



1 + exp(-2/3a;) 



(1) 



the stock state Si(t) is reinforcing or persistent so that 
Si(t + 1) = Si(t), while with probability 1 — Q(uj(t)) it 
reverses or is anti-persistent so that Si(t + 1) = — Sj(t). 

Therefore it is a basic assumption of this model that 
there are always persistent and anti-persistent stocks and 
we can imagine how this might occur as follows: Every 
time t, the speculators choose a pair of stocks a and [3 
from, say, the excess demand stocks so that s a (t) = 1 and 
sp(t) = 1 and compare them by comparing their confi- 
dences V a (t) and Vp(t). Suppose V a (t) > Vp(t), therefore 
the speculators decide to demand a more and to supply 
(3 instead. Therefore s a (t + 1) = 1 and sp(t + 1) = —I. 
Conversely suppose the speculators are comparing stocks 
a and j3 which are in excess supply with s a (t) = —1 
and Sjs(t) = —I, where V a (t) > Vp{t), then speculators 
decide to continue to supply a while demanding (3 so 
that s a (t + 1) = — I and sp(t + 1) = 1. If we imagine 
this process of comparing stocks happening continuously 
we can imagine each share generally interacting with the 
market confidence V(t) through the relative confidence 
Ui(t) = Vi(t) — V(t) and Eqn|y with the above Sj(t) up- 
date dynamics becomes a plausible description. (In fact 
as an alternative model, to be reported separately, we can 
say high valued stocks a are persistent s a (t + 1) — s a (t) 
but instead of low- valued stocks (3 being anti-persistent we 
could imagine that they are completely independent of the 
previous state such that s@(t + 1) = ±1 at random.) (We 
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Fig. 1. Single realization of the market index price time series 
p{t) for c = 0.001, P = 80 and N = 200. 



should also point out that if the idea of anti-persistent 
reversing demand/supply stock states seems strange so 
should the idea of persistent high volatility itself!) 

The relative confidence Ui(t) measures to what extent 
speculators have faith in the persistence of stock i's state 
Si(t) compared to their belief in the persistence of the 
whole market V(t). The 'inverse temperature' parameter 
(3 in Eq.0 measures to what extent speculator expectation 
is self-fullfilling, that is the extent to which homogenous 
'herding' occurs. When (3 = oo the system is determinis- 
tic and Eq.fj] reduces to a step function and the state of 
stocks i, Si[t), with V-i{t) > V(t) will be persistent with 
probability 1, while stocks i with Vi(t) < V(t) will be 
anti-persistent with probability 1, i.e. complete homoge- 
nous herding. In this case the speculators behave with 'one 
mind'. With (3 = all states Sj(i+1) are chosen randomly 
independent of Si(t) and speculators are therefore com- 
pletely independent in their viewpoints. (We can also say 
equivalently that (3 measures the extent to which traders 
know stock i's confidence Vi{t).) It is a basic assumption 
of this model that a single stock cannot be considered on 
it's own, independently of the rest of market, speculators 
are always dynamically choosing. This is reminiscent of 
an ecosystem where a single organism cannot be consid- 
ered independently of the other coevolving organisms in 
the system. 

In the version of this model presented here, the dy- 
namics of the confidences Vi (t) themselves are treated as 
internally defined behaviour of the model. If a stock i has 
persistent excess demand/supply at time t then we sug- 
gest speculators decrease their confidence V,(t) in it by a 
small amount c, (c > 0). If on the other hand the excess de- 
mand/supply of stock i is observed to reverse we say spec- 
ulator confidence is coupled to the excess demand/supply 
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state of the market G(t) such that the change in confi- 
dence, AVi(t) = Vi(t + 1) - Vi(t), is given by, 



AVi(t) = 



— c, if 
Si(t)G(t), if 



s t (t + l) = Si(t) 
s t {t + l) = - Si (t) 



(2) 



Therefore since Si(t)G(t) — —Asi(t)G(t)/2, specula- 
tor confidence Vi(t) in stock i decreases when stock z's 
state Si(t) moves into the overall market majority state as 
measured by G(t). This means the probability of a stocks 
state reversing increases when it moves into the majority 
and decreases when it moves into the minority. This rule 
has some similarity to the Minority Gamefl"^]. We believe, 
however, that it is more appropriate for financial market 
modeling to increase an agent's fitness when the agent 
moves into the minority, rather than when it is in the 
minority. As is well known, the way to make speculative 
money in markets is to stay in the majority but switch 
state into the minority just before everybody else does. 
(See alsoQ.) Of course in the Minority Game the agents 
are traders not stocks, but here we assert that a similar 
effect occurs between individual stocks and the market it- 
self because of the effect of the coupling of the individual 
stocks to the overall market. This coupling is due to the 
collective choosing action of the background traders. 

Our stock confidences Vi (f ) and the corresponding prob- 
abilities Qi(t) are defined in cumulative terms, through 
continuous changes, like cumulative fitnesses in ecody- 
namics. We believe this is in keeping with the way people 
think, traders do not simply forget their previous evalu- 
ations of the stocks, constructing probabilities anew, but 
rather dynamically update their perceptions, in a simple 
way depending on the current behaviour of the stock. 

While the parameter (3 measures the amount of homo- 
geneity of speculator opinion in the market, c measures 
the rate by which the speculators in general lose faith in 
a current price trend Sj(i + 1) = Si(t). That is the rate by 
which they lose faith in a stock's bullishness or bearish- 
ness. 

The dynamic is synchronous, all Vi(t) and Sj(i) are up- 
dated at the same V(t) and G(t) according to Q(ui(t)). For 
initial conditions the spins Sj (0) are chosen randomly and 
the confidences Vi (0) randomly and uniformly on [—1,1]. 



3 Numerical Results. The Behaviour of the 
Volatility 

In our paper ||] we describe the behaviour of the model in 
detail as we vary the parameters c and (3. As explained in 
that paper we believe real financial markets are described 
by c « (c > 0) and (3 « 80. We show in that paper 
that for those parameters values the price returns distri- 
butions agree very well with the empirical observations 
mentioned in the Introduction. A typical index price time 
series for c = 0.001, = 80 and N — 200 is shown in 
Fig.|l|. This is just defined from the cumulative changes as 
p(t) = Ylj—i Ap(J), Qualitatively it looks very reminis- 
cent of a real time series, with periods of time which look 
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Fig. 2. Distribution of absolute price changes |/lp(t)At| for 
At = 9, c = 0.001, (3 = 80 and N = 200. The distribution is 
calculated from one single realization of the price index time 
series. The lower panel shows the distribution in linear-linear 
axes, the upper panel in log-log axes. 



like Gaussian random walks and other times which have 
larger fluctuations. One 'crash' is visible, as is a period of 
'sustained growth'. 

Here we study the behaviour of the volatility for the 
same parameter values and compare it with real observa- 
tions mentioned in the Introduction. 

Volatility is usually measured from financial market 
time series by taking a certain time window and calcu- 
lating the standard deviation of the price fluctuations in 
that window. When considering time series of length sev- 
eral years the window is often of the order of a few weeks. 
Here, as a proxy for the volatility, we use the distribu- 
tion of the absolute changes in price in a similar way 
to[^2). This is shown in Fig.|| for the absolute values of 
Ap(t) At=9 = X^tt Ap(j). We choose to show the At = 9 
distribution, however other At distributions are similar. 
Shown in the bottom panel is the distribution in linear- 
linear axes. The distribution is similar in appearance to 
the log-normal distribution, although it has a long tail. 
The distribution is shown in log-log axes in the upper 
panel of Fig.||. The tail is a straight line meaning the ex- 
treme absolute changes have a power-law distribution. The 
slope of the tail is approximately 4, in very good agree- 
ment with empirical results in |1^ | , explained in the Intro- 
duction. 

Next we study the dynamical behaviour of the volatilty. 
In |j) we show that the price changes time series from this 
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Fig. 3. Diffusion of the volatility for a single realization of the 
price index time series with c = 0.001, f3 — 80 and N = 200. 
The variance a 2 (At) is plotted against At in log-log axes. The 
upper curve is the results from the time series, the lower curve 
is from the scrambled time series. 



model show intermittent volatility clustering. In that pa- 
per we give some explanation for this behaviour. For our 
chosen parameters c = 0.001 and j3 = 80 the time se- 
ries shows weak superdiffusion, the characteristic of per- 
sistence. Shown in Fig in log-log axes is the variance 
a 2 (At) of the price changes Ap(t)^t time series, plotted 
against At. From the relationship, 



a 2 (At) = (At) 



2/i. 



(3) 



we obtain 2/i w 1.5 for At rs 1 ~ 300, while for At > 
300 we obtain 2/i w 1.0. This implies superdiffusion at 
shorter time scales but with reversion to normal diffusion 
for longer time scales. This is in very close agreement with 
studies in |l2| , explained in the Introduction. Also shown 
in Fig.|| is the same analysis for the scrambled time series, 
performed as a check. As expected the slope reverts to 
that for normal diffusion. 



4 Discussion 

We have described a multi-share financial market model 
and numerically studied the behaviour of the volatility of 
the market index price changes Ap(t) for the parameter 
values c = 0.001 and (5 — 80. Our results agree very well 
with empirical studies. In our paper j^] we show these re- 
sults are robust to changing the value of c providing it 
remains 'small' or equivalcntly the system size N is large. 
Since c measures the rate at which speculators become 
nervous of continuing bullish or bearish trends, c — > 0, 
(c > 0), means markets find a state where speculators try 
to collectively 'stand their ground' the longest. (3 ss 80 in 



fact corresponds to a phase transition between a regime 
where price fluctuations obey Gaussian statistics and a 
regime where they obey Levy distribution [|| statistics. 
In 0] we give qualitative reasons why markets should be 
found in this phase transition region. 
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